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Cloning and expression profiling of heat shock protein 
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Abstract; [ Aim] The sHSPs are a diverse of molecular chaperones with molecular weights ranging from 12 to 
43 kDa in different species. These proteins have been shown to be crucial components of the diapause of some 
species, but need to be further elucidated. The study aims to clone and characterize the sHSP gene of Delia 
antiqua, of which orthologs were earlier reported to be diapause-related, and to investigate the expression 
patterns of the gene during both winter and summer diapauses of the species. [Methods] A sHSP gene was 
cloned from D. antiqua using RACE-PCR, and characterized through homology, domain and phylogenetic 
analyses. The existence of the intron was investigated with genome sequencing around the gene. The expression 
level of this gene during both winter diapause (WD) and summer diapause (SD) was detected real-time qPCR, 
and the expression difference was compared to reveal the association with diapause development. [Results] The 
cDNA of sHSP cloned from D. antiqua, which is named DaHSP23 (GenBank accession no. HQ392521.1), is 
904 bp in length, encoding 186 amino acids with a molecular mass of 20.9 kDa and a theoretical isoelectric 
point of 6.42. DaHSP23 shares more than 66% amino acid identity with the known HSPs of other dipteran 
insects and is homologous with the known HSP23 genes of other dipteran insects. Genome sequencing around 
the gene revealed that it is intronless. The expression of DaHSP23 was up-regulated in both SD and WD pupae, 
with higher expression level detected in the late maintenance stage, and down-regulated when those diapauses 
were terminated. 【 Conclusion] The patterns of HSP23 regulation are substantially different between different 
developmental stages of both SD and WD, and its regulation in diapause might be species-specific. DaHSP23 
might play an essential role in these two types of diapauses in D. antiqua. 
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1 INTRODUCTION 


Heat shock proteins (HSPs) are conserved 
proteins involved in multiple cellular processes, 
including protein folding, targeting, and translocation 
across membranes ( Neupert, 1997; Hartl and Hayer- 
Hartl, 2002 ). 
molecular weights, HSPs can be mainly classified into 
three main families: HSP90 (85 - 90 kDa), HSP70 
(68 -73 kDa) and the small HSP (sHSP) (12 -43 
kDa) (Denlinger et al., 2002). 

sHSPs contain an aq-crystalline domain with 
molecular weights of 12 -43 kDa, depending on the 
variable N- and C-terminal extensions ( MacRae, 
2000; Taylor and Benjamin, 2005). Studies in 


According to the homology and 
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different experimental systems have revealed a 
variety of functions for the sHSP under stress 
conditions, including basic chaperoning activity 
(Ehrnsperger et al., 1997; Haslbeck et al., 1999) , 
cytoskeleton protection ( Lavoie et al., 1993) and 
modulation of the apoptotic process ( Bruey et al., 
2000). Despite the classical definition of heat shock 
proteins as polypeptides induced by stress, cell- 
specific expression of sHSPs in the absence of stress 
has been reported during the development of a wide 
range of organisms such as Caenorhabditis elegans 
(Ding et al., 2000), Drosophila melanogaster 
(Glaser et al., 1986; Marin et al.,1993; Michaud et 
al.,1997) , Xenopus laevis (Lang et al., 1999) , Mus 


musculus ( Gernold et al., 1993; Tanguay et al., 
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1993; Armstrong et al., 2001), and human being 

( Jantschitsch et al., 1998). 
Diapause was developed by 

arthropods to survive from unfavorable seasons. Delia 


insects and 


antiqua can be induced into both summer diapause 
(SD) and winter diapause ( WD ) after head 
evagination in pupal stage, and easy to be reared in 
laboratory; therefore, it is regarded as a good model 
for diapause study (Ishikawa et al., 1987, 2000). 
The diapause program was suggested to associate with 
the expression of diapause-specific 
( Denlinger, 2002 ). 


genes 
The heat shock proteins, 


including sHSP, HSP60, HSP70 and HSP90 
(Macrae, 2010), are of great importance in 
understanding the molecular mechanisms of 
diapause. 


In response to temperature and stresses, HSP90 
and HSP70 were expressed in the summer and winter 
diapause pupae of D. antiqua (Chen et al., 2005a, 
2005b, 2006). In the pupal diapause of the flesh fly 
Sarcophaga crassipalpis, 23 kDa heat shock protein 
(HSP23 ) gene transcripts are expressed at very high 
levels as part of the normal diapause program, 
indicating that HSP23 plays an essential role during 
this over-wintering developmental arrest ( Yocum et 
al., 1998 ). 


sericata, the level of HSP23 transcripts was found to 


However, in the blow fly Lucilia 
be consistently low, and no difference in its 
expression level was observed between normal 
condition and diapause , indicating that expression of 
HSP23 is not regulated in response to diapause 
(Tachibana et al., 2005). This inconsistency needs 
to be elucidated with additional species and different 
types of diapause. To our knowledge, no study has 
as yet been carried out to explore the expression of 
HSP23 during insect summer diapause. It will be 
interesting to compare the expression patterns 
between SD and WD. 

In the present study, we cloned the full-length 
cDNA of HSP23 gene ( DaHSP23) from the onion 
maggot, D. antiqua, compared the expression 
patterns of DaHSP23 among summer-, winter- and 
non-diapausing pupae using real-time quantitative 
RT-PCR, and discussed the HSP23 regulation in 


relation to diapause. 


2 MATERIALS AND METHODS 


2.1 Insect rearing and diapause induction 

The non-diapausing (ND) D. antiqua was 
reared on an artificial diet at 20°C with a 16L: 8D 
photoperiod and 50% - 70% relative humidity as 
described by Chen et al. (2010). To induce SD, 
larvae were maintained at 25°C with the 16L: 8D 


photoperiod. Newly formed pupae were kept under the 
same conditions as larvae until 10 days after pupation 
(D10) and then transferred to 16°C and 16L: 8D 
photoperiod to trigger diapause termination (Chen et 
al., 2010). To induce WD, larvae were reared at 
15°C. with the 12L: 12D photoperiod (Ishikawa et al., 
2000), and pupae were kept under this condition 
throughout. Pupae collected at different developmental 
stages of ND, SD and WD were frozen in liquid 
nitrogen and stored at —80°C for RNA extraction. 
2.2 RNA extraction, cDNA synthesis, and 
PCR amplification 

Total RNA was isolated with the RNeasy Mini 
Kit ( QIAGEN ) supplier ” s 
instructions. The first-strand cDNA was synthesized 
from total RNA using oligo (dT) primer (TaKaRa) 
and Superscript™ reverse transcriptase ( Invitrogen ) 


according to the 


according to the supplier’ s instructions. A 284 bp 
cDNA fragment of HSP23 was amplified by PCR from 
cDNA using degenerate primers ( Table 1), which 
were designed based on consensus sequences ( number 
from 202 to 485, Fig. 1.) from several dipterans. 
PCR conditions were as follows; 1 cycle at 94°C for 3 
min; 35 cycles at 94°C for 30 s, 55°C for 30 s and 
72°C. for 60 s; and 1 cycle at 72°C for 10 min. The 
expected DNA band was excised from the agarose gel, 
purified using a DNA gel extraction kit ( Bioteke, 
a pMD-19T Vector 
Plasmids were isolated from positive 


Beijing ) and cloned into 
( TaKaRa ). 
clones and sequenced by Shanghai Songon Company. 
2.3 Amplification of the full-length HSP23 cDNA 

To obtain the full-length HSP23 cDNA, we 
performed 5’ and 3' rapid amplification of cDNA 
ends (5’ RACE and 3’ RACE, respectively) using 
SMART™ RACE kit (Clontech) according to the 
manufacturer’ s instructions. For 5' RACE, gene- 
specific primers Gene Races 5'-R and Gene Races 
5n’-R were designed based on the 284 bp cDNA 
fragment. The amplification was performed with 1 
uL of 5’-ready-cDNA using Universal Primer Mix 
(UPM, Clontech) and Gene Races 5'-R, and then 
the nested PCR was carried out with Gene Races 
5n’-R and Nested Universal Primer A ( NUP, 
Clontech). The PCR conditions were as follows; 1 
cycle at 94°C for 3 min; 35 cycles at 94°C for 30 s, 
at 60°C for 30 s, and at 72°C for 60 s; and 1 cycle 
at 72°C for 10 min. For the 3’ RACE, the cDNA was 
synthesized according to the manufacturer’ s protocol 
[ TaKaRa RNA PCR Kit (AMV) Ver. 3.0]. Gene- 
specific primers Gene Races 3’-F and Gene Races 
3n'-F were designed also based on the 284 bp cDNA 
fragment. The amplification was performed with 
1 uL of 3’'-ready-cDNA using M13 Primer M4 (M4, 
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Primer name 


Table 1 Primers used in this study 


Primer sequences (5' -3’) Use of primers 





Hsp23F 
Hsp23R 
Gene Races 3'-F 
Gene Races 3n'-F 
Gene Races 5'-R 
Gene Races 5n'-R 


TTYCARGTNTGYATGGAYGT 
GGNCSNACYTGYTGDAYTG 
CGAGGGTAAACATGAGGAACGGGAA 
ACCCGACACTTTGTACGTCGCTATGC 
GCATAGCGACGTACAAAGTGTCGGG 
TTCCCGTTCCTCATGTTTACCCTCG 


GAA GAA AAC AAA ACC AAAGCA AAT CAA ATC 
AAC AAA GTC AAA TTA ATA TTT TTA TIT AAA AAC TTT AAATAA TTG TTA AAA 


Hsp23-CDSF ATGGCAAACTTACCACTTCTTTTGA 
Hsp23-CDSR TTACAATTTCGATGGAGGCTGTTCC 
Hsp23QF ACACTTTGTACGTCGCTATGCTT 
Hsp23QR ACGTTCGTTGCTGGATTTATCT 
-132 
-102 
-51 


1 
1 


52 
18 


103 
35 


154 
52 


205 
69 


256 
86 


307 
103 


358 
120 


409 
137 


460 
154 


S511 
171 


562 
613 
664 
715 
766 


AGT GTT GAA ATT TTA AAT ATT ATT AAA TAA TAA AAC AGT AAA TCT ATA AAA 
ATGGCA AAC TTA CCA CTT CTT TTG AGT TTG GCT AAC GAT CTC AAT CGC CTT 
M A N L P L L L S L-A N DBD L N R L 


TCC TTG GCT GCT TCG CCG TTT TAT GAATCA CCA ATT CAC TAT TTG CGC CAT 
S L A A S P F Y E Ss FP I H Y L R 4H 


CCT GGC TAT TTG GCC TTG CTG AAC GCT GAT CAA CCA CAACAA CTA CGAAAA 
Bo Ge eo A a. ENT SAD ae BY SG) s@? ay eRe ie 


CAT GAT AAGGAATCG TCT GGT CCA CTA GCC ACC GTT GGC AAG GAT GGT TTC 
H D K E S S G P L ATE o SG o ReDD G F 


CAG GTA TGC ATG GAT GTT CAA CAA TTC AAG CCC AGT GAG TTG AAT GTC AAA 
人 


GTC ATC GAC AAT TGT GTA ATT GTC GAG GGT AAA CAT GAGGAA CGGGAA GAT 
Vv I DN C Vv I Vv E G K H E E R E D 


GAT CAT GGC TTT ATT ACC CGA CAC TTT GTA CGT CGC TAT GCT TTA CCC AAA 
D H G F I T R ë H F Y È R Y A L P K 


GGT TAT GAT TCC AAT AAA GTT CAA TCC ACT TTA TCC TCG GAT GGT GTT TTA 
C ae na S N K V Q S T L S S D G V L 


ACA GTT AGC GTA CCC AAA CCT CAA ATC GAA GAT AAA TCC AGC AACGAA CGT 
T Vv S ki P K P Q I E D K S S N E R 


CAA ATT CAA ATT CAA CAA GTT GGT CCA GCA CAT TTG AAT GTC AAAGAA AAT 
CT e AT: SGN Cor Gaye Ge SR A ET Gi A See E N 


CCT AAA GAA GAAAAG GAA AGT GAT AAG GAA CAG CCT CCA TCG AAA TTG TAA 
P K E E K E S D K E Q P P S K ë L x 


AAGTCT AAC ACT TAT TGGCAT TTG TTT GAT TAT TTG AGAGTA AAT AGT CTG 
TTT AAT TTT AAGTAA AAA TTT TAA TTA TAC TTG TTA ACC AAA ATT TTA TAC 
GAC TAT TTA TAA AAA TTA AAAAGA ACT GAA TTT TGT TTT AAA AATAAATAA 
AAT TTA TAA GTA CCA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAAAAAAAA 
AAAAAA A 


Fig. 1 cDNA and deduced amino acid sequences of HSP23 from Delia antiqua ( DaHSP23 ) 


The start codon ATG and stop codon TAA are indicated in bold. A putative polyadenylation site (AATAAA) is underlined. The numbers on the left 


indicate the positions of nucleotides and amino acids in the sequences. 


TaKaRa) and HSPF2, and then the nested PCR was 
carried out with HSPf22 and M13 Primer M4 (M4, 2.4 Cloning of HSP23 gene from D. antiqua 


TaKaRa). The PCR conditions were as follows; 1 genomic DNA 


cycle at 94°C for 3 min; 35 cycles at 94°C for 30 s, at 


Amplification of DaHsp23 fragment 
Amplification of DaHsp23 3’ end 
Amplification of DaHsp23 5' end 

Amplification of DaHsp23 genomic DNA sequence 


DaHsp23 real-time qPCR 


(http://www. mbio. ncsu. edu/ BioEdit/bioedit. html) . 


59°C for 30 s, and at 72°C for 30 s; and 1 cycle at 
72. for 10 min. 

The obtained 5’- and 3'-RACE PCR products 
were cloned, sequenced and assembled into the full- 
length HSP23 cDNA sequence using BioEdit software 


Genomic DNA was isolated from non-diapausing 
pupae of D. antiqua using the Universal Genomic 
DNA Extraction Kit Ver. 3.0 (TaKaRa, Japan) 
according to the manufacturer’ s instructions. The 
concentration and quality of DNA were verified by 
spectrophotometer and electrophoresis on 1. 0% 
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agarose gel. In order to determine whether the 
DaHSP23 gene contains introns in the coding region, 
the genomic DNA sequence for DaHSP23 was 
amplified from genomic DNA (100 ng). The 
amplification was performed using a pair of specific 
primers Hsp23-CDSF and Hsp23-CDSR, and the 
amplification conditions were as follows; 1 cycle at 
94°C for 3 min; 35 cycles at 94°C for 30 s, at 59°C 
for 45 s, at 72°C for 1 min; and 1 cycle at 72°C for 
10 min. The genomic DNA sequence was then 
cloned and sequenced for DaHSP23 intron analysis. 
2.5 Sequence and phylogenetic analysis 

The nucleotide and protein sequences were 
edited, and the theoretical molecular weight and 
isoelectric point of sHSPs were predicted using 
Biokdit software ( http://www. mbio. nesu. edu/ 
BioEdit/bioedit. html ) . 


sequence of the DaHSP23 cDNA was used as query 


The deduced amino acid 


to perform a global search in GenBank non-redundant 
protein database by BlastP program at the NCBI 
website (http://www. ncbi. nlm. nih. gov’BLAST) , 
and then the dipteran homologues of sHSPs were 
selected and analyzed. Multi-sequence alignment was 
generated using Clustal X and DNAStar 5.0 software. 
The functional motif and sites were deduced in 
reference to homologous sHSP of D. melanogaster. 
The phylogenetic analysis using neighbor-joining 
method with p-distance model was performed with 
MEGA 4.0. The homologous sHSP of 
Caenorhabditis elegans ( NP_001023957. 1) was 
used as the outgroup, and a total of 1 000 bootstrap 
replications were used to test the topology. 
2.6 Real-time qPCR 

Total RNA was isolated from 30 mg of pupae 
using a RNeasy Midi/Maxi kit (Qiagen) and treated 
with RNase-Free DNase I ( Qiagen). cDNA was 
subjected to reverse transcription from total RNA 
using random primers with a RNA PCR kit 
(TaKaRa) and purified with a SUPREC -02 kit 
(TaKaRa). The HSP23 fragment (133 bp) was 
amplified using gene-specific primers Hsp23QF and 
Hsp23QR. The 18S rRNA gene was used as a 
reference. Primers for 185 rRNA were the same as 
those used by Chen et al. (2005a) , which amplifies 
a 333 bp fragment. One cDNA sample, diluted to 1, 
571,57, 57, 57 and 5~°, was employed as an 
internal standard. Real-time qPCR was performed in 
25 uL reactions containing 1 L of template cDNA 
or the standard, 2x SYBR Green PCR premix 
(TaKaRa), and 0. 3 mmol/L of each primer. 
Thermal cycling conditions were as follows; 1 cycle 
at 94°C for 15 s, 40 cycles at 94°C for 10 s, 60°C 


for 30 s. Melting curve analysis was performed from 


60 to 94°C. The relative molar levels of HSP23 and 
18S rRNA transcripts were calculated based on 
crossing point analysis, using standard curves 
generated from the cDNA standards. HSP23 mRNA 
levels were normalized with those of 18S rRNA in the 
same samples quantified in the same manner, and 
the final relative mRNA levels of HSP23 were 
averages of three replicates. The fluorescent real- 
time PCR assay was carried out in a Bio-Rad iCycler 
iQ5 ( Bio-Rad, USA ). Statistical analysis was 
performed using SPSS 13.0 for Windows. Data were 
expressed as mean + SE and subjected to a one-way 
ANOVA, followed by Duncan’ s test and a Student’ s 
test. A P <0.05 was considered to be significant. 


3 RESULTS 


3.1 Cloning and sequence analysis of DaHSP23 

A fragment of 284 bp was amplified using the 
degenerate primers, and its nucleotide sequence is 
significantly similar to those of other known HSP23s 
in GenBank (E value, 8e~° —3e~’). The 94 a. a. 
sequence encoded by the cDNA also has high 
similarity to known HSP23 protein sequences ( E 
value, 5e” - S5e™). the cDNA 
fragment was determined to be the partial sequence 
of HSP23 gene. 
obtained, the gene-specific primers were synthesized 
for 3’-end and 5'-end RACE. And then, a 500 bp 
fragment was amplified by 3’ RACE and a 480 bp 
fragment by 5’ RACE. The full-length cDNA 
sequence of DaHSP23, which has been deposited in 
GenBank under accession number HQ392521. 1 and 
ADX36150.1, was obtained by alignment analysis of 
the above three fragments (Fig. 1). 


Therefore , 


Based on the cDNA sequence 


For DaHSP23 intron analysis, the genomic 
DNA sequence for DaHSP23 was amplified by 
specific primers (CDS23F, CDS23R), and then 
was compared to cDNA sequence for DaHSP23. The 
result showed that these two sequences are identical , 
indicating that there is no intron in the coding region 
of DaHSP23 gene. 

The full-length DaHSP23 cDNA, 904 bp in 
length, encoding a 186 a. a. peptide with a calculated 
molecular weight of 20. 9 kDa and a theoretical 
isolectric point of 6. 42, contains a 5’-terminal 
untranslated region (UTR) of 132 bp, a 3’-terminal 
UTR of 211 bp and an open reading frame (ORF) of 
561 bp. A polyadenylation signal sequence AATAAA 
is located at nucleotides 710 - 715 (Fig. 1). The 
alignment of DaHSP23 with other five homologous 
dipteran sHSPs ( DmHSP20. 6, DbHSP20. 6, 
CcHSP18. 9, GmHSP23. 1 and ScHSP23. 0) is 


shown in Fig. 2. The conserved q-crystallin domain 


406 昆虫 学 报 Acta Entomologica Sinica 59 卷 





is predicted to consist of a. a. 65 to 142, which is 
The InpA 
posttranslational modification , 


characteristic of sHSPs. region for 
protein turnover and 


chaperones is predicted to consist of a. a. 12 to 158, 










DaHSP20.9 MANLPLLLSLANDLENRLSLAASPFYES---- 
CcHSP18.9 MANLPLILSLADDLSRLS-—--PLYE----- 
DbHSP20.6 MANLPLLLSLADDLNRLS—-MVPFYDT——-- 
DmHSP20.6 MANIPLLLSLADDLGRMS--MVPFYE----- 
GmHSP23. 1 
ScHSP23. 0 

| ag | 90 100 
DaHSP20.9 -------- GPLATYGRDGFQ 
CeHSP18. 9 --------— PIASVGKDGF 
DbHSP20.6 ----—- SIGPVSKIGKDGFQ 
DmHSP20.6 ASSG--SSGAVSKIGKDGFQ 
GmHSP23.1 RRGEKRGEDRVSAIGKDGFO 
ScHSP23.0 RRGGEKKEGRVSPYGKDGFQ 

| un 1 | 170 
DaHSP20.9 GYDSNKVQSTLSSDE 
CcHSP18.9 GFEADKVYSTLSSD 
DbHSP20.6 GYEADKYASTLSSD 
DmHSP20.6 GYEADKVASTLSSD| 
GmHSP23.1 GYDADKVVSTLSSD 
ScHSP23.0 GYDADKYVSTLSSD 









PKPSELNVRY I DNCVI VEGKHEEREDDHGF IT RHF 
FKP SELSVKVVDDHI VVEGKHEEREDDHGY ISRHFY y 
FKP SELNAKVVDNS I VVEGKHEEREDDHGF ITRHFVR 
SHFKPSELY VK VQDNSVL VEGNHEEREDDHGFITRHFYRRY 
OF MM Wie o RY 


based on D. melanogaster ( AAN11962. 1). The 
putative dimer interface sites are predicted to be at a. 
a. 71 (C), a.a. 73 (D), a.a. 75 -76 (QQ), a.a 

113 (R), and a.a. 134-135 (GV) (Fig. 2). 


—---PYYWQRQNRNP---YLALAAVPLDQQL ROLEKQVG 
——--PYYCQRQ-RNP---YLAL VG—PMEQOLRQLEKQVG 
MSSLPLLLSLIDELDREANH—-PSYHGNDF GLGLPPYL IHRHPYREPANNMIGYSLPFGI ISRLNERQAG 
MAGLPLLLSLIDELDRESN---PSYYGNDFGLGLSPYL IHRQPQREPAHNL VGYSLPLSLLSRLNEHQVA 


110 120 130 140 
clear eek wee ace eset ca 
JARTA LPK 
RSYALPK 
RYALPD 







[TVS¥PRP-O1EDKSSNERO1Q1QQVGPAHUNVKENPKEEKESDKEQPPSKL—-—— 
eh wena NeRv vA QT ERAHUNVKENPKD--AANEEKPKA 一 一 -一 


vs YPR QTEpRA Belk VR ee ok vcDAANGGGENAEARA 
LTVS ¥PKP—QEEEKS-NERV 1 Q1QQT GPAHLNVKONT QEK VGDGANGEPQ-EEKKA 


Fig. 2 Alignment of DaHSP23 (DaHSP20.9 in the figure) with other five homologous dipteran sHSPs 


The conserved q-crystallin domain is underlined, the InpA region is double underlined, and the putative dimer interface sites are boxed. The amino 


acids with 100% and over 70% identity are shaded with dark and gray color, respectively. A “ - ” corresponds to a gap and the residue column is each 


assigned a number above. The accession number of each sequence is same as in Fig. 3. 


3.2 Homology and phylogenetic analysis of 
DaHSP23 


The global search with BlastP in GenBank non- 
redundant database using DaHSP23 as query 
retrieved 18 homologues with E-value larger than 
4e ”and 100 with E-value larger than le”. Of the 
100 homologues, 28 belong to Diptera, and 10 
representative sequences (Table 2) with full length 
were finally selected to perform further sequence and 
phylogenetic analysis with DaHSP23. The DaHSP23 
showed the highest identity ( 68. 2%) with 
DvSHP20.7 (Dr. virilis) and DbHSP20. 6 ( Dr. 
buzzatii), followed by CeHSP18. 9 (65. 1%, 
GenBank accession no. ACG58884. 1), DmHSP20. 6 
(64.9% , GenBank accession no. NP_523999. 1) 
and DpHSP20. 4 (64.4% , GenBank accession no. 
XP_002024748.1) (Table 2). 

The phylogenetic tree of DaHSP23 and the other 
10 dipteran homologous sHSPs was constructed, with 
the homologous sHSP of C. elegans as the outgroup 
(Fig. 3). Six sHSPs from six species of Drosophila , 


two sHSPs of Ceratitis capitata, and two sHSPs from 
Glossina 


morsitans morsitans and Sarcophaga 
separately 
monophyletic clades, all with 100% bootstrap support. 
3.3 Expression of DaHSP23 during summer 
and winter diapause 

The relative mRNA levels of DaHSP23, 
normalized to 18S rRNA levels, were very low in 
both pre-diapause (0. 85 and 1. 46 for SDO. 5 and 
SD1, respectively ) and post-diapause (2. 19 for 
SD17) stages (Fig. 4). At the onset of SD (SD5), 
the level increased to 3.29, and then rapidly peaked 
(4. 50, SD10). Thereafter, the level gradually 
decreased from 6.04 (SD11) to 2.19 (SD17), with 
a slope rate of -0. 64/d (Fig. 4). The relative 
mRNA levels of DaHSP23 , normalized to 18S rRNA 
levels, were relatively low in the pre-diapause stages 
[0.97(WD3), 1.06( WD4) |. During the diapause , 
the level rapidly increased from 2.98 ( WD29) to 
8.01 (WD54) , and then gradually decreased to 2. 33 


(WD103), and then to 2.22 (WD114) , with a slope 


crassipalpis were clustered into three 
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Table 2 The dipteran sHSPs homologous to DaHSP23 and used for sequence and phylogenetic analysis , 
and their identity and similarity compared with DaHSP23 
dikes Sadie GenBank Amino acid sequence Amino acid sequence 
accession number identity (% ) similarity ( % ) 

DmHSP20. 6 Drosophila melanogaster NP_523999. 1 64.9 74.7 
DpHSP20. 4 D. persimilis XP_002024748. 1 64.4 74.7 
DwHSP20. 8 D. willistoni XP_002062617. 1 61.5 76.6 
DvSHP20. 7 D. virilis XP_002048201. 1 68.2 79.7 
DbHSP20. 6 D. buzzatii ABX80640. 1 68.2 79.2 
DgHSP20. 4 D. grimshawi XP_001985085. 1 66.0 77.3 
CeHSP18.9 Ceratitis capitata ACG58884. 1 65.1 79.6 
CcHSP18.9 C. capitata ACG58883. 1 62.9 77.4 
GmHSP23. 1 Glossina morsitans morsitans ADD18976. 1 53.1 66.4 
ScHSP23. 0 Sarcophaga crassipalpis AAC63387. 1 57.0 69.1 
Outgroup Caenorhabditis elegans NP_001023957. 1 30.3 45.7 


sHSP sequence names are suffixed with deduced molecular weights in kDa. 
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DmHSP20.6 | NP_523999.1 


DpHSP20.4| XP_002024748.1 


DwHSP20.8 | XP_002062617.1 


DvSHP20.7 | XP_002048201.1 
DbHSP20.6 | ABX80640.1 
DgHSP20.4 | XP_001985085.1 
DaHSP20.9 | ADX36150.1 $ 
CcHSP18.9 | ACG58883.1 


CeHSP18.9| ACG58884.1 


GmHSP23.1| ADD18976.1 


ScHSP23.0 | AAC63387.1 


Caenorhabditis elegans | NP_001023957.1 


Fig. 3 The neighbor-joining (NJ) tree of eleven dipteran sHSPs 


The homologous sHSP of Caenorhabditis elegans (GenBank accession number; NP_001023957.1) was used as the outgroup, and percentage bootstrap 
virilis; Db: D. buzzatii; Dg: D. 


values are indicated on each node. Dm; Drosophila melanogaster; Dp: D. persimilis; Dw: D. willistoni; Dv; D. 
grimshawi; Ce; Ceratitis capitate; Cc; C. capitate; Gm; Glossina morsitans morsitans; Sc; Sarcophaga crassipalpis. 


rate of —0.09/day (Fig. 4). In brief, the relative 
mRNA levels of DaHSP23 in the diapause stages 
(WD29, WD54, and WD79) were much higher than 
those in pre-diapause ( WD3, WD4 ) and post- 
diapause stages ( WD103, WD114) (P < 0. 05) 
(Fig. 4). These results indicated that the expression 
of DaHSP23 was almost constitutive under unstressed 


conditions but was rapidly up-regulated during winter 
diapause. The relative mRNA levels of DaHSP23 in 
SD and WD were similarly high compared to those in 
ND (between 0.94 at ND1 and 1.46 at ND9) (Fig. 
4). The relative mRNA levels of DaHSP23 in pre- 
diapause and post-diapause of SD and WD were 
similarly low compared to those in ND (Fig. 4). 
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Fig. 4 Real-time qPCR analysis of DaHSP23 mRNA levels in 

non-diapausing (ND) (A), summer diapausing (SD) (B) 

and winter diapausing ( WD) (C) pupae of Delia antiqua 
Relative mRNA levels were normalized to the internal standard 18S 
tRNA, with the bars indicating the mean + SE of three replicates and 
different letters between samples demonstrating the significant transcript 
difference ( P < 0. 05). Data were expressed as mean + SE and 
subjected to a one-way ANOVA, followed by Duncan’ s test and a 
Student’ s test. ND1, SD0. 5 and WD3 represent day 1, day 0.5 and 
day 3 after pupation of ND, SD and WD, respectively. 


4 DISCUSSION 


The amino acid sequence of DaHSP23 contains 
an alpha-crystallin domain ( ACD) and a hallmark, 
which is a characteristic of the sHSP family ( De Jong 
et al., 1998; McHaourab et al., 2009 ). The 
sequence exhibits high conservation with other 
dipteran HSP23s, and shares the highest similarity 
with those reported to be related to diapause. 

In the HSP70 family, introns are generally 
found only in constitutively expressed genes but not 
in the inducible genes ( Ming et al., 2010). In sHSP 


family, most members are intronless ( Caspers et al., 


1995). It was reported that the lack of introns in 
inducible HSP genes could help to circumvent a 
block in RNA splicing and enable preferential 
expression of the heat shock proteins during periods 
of stress ( Yost and Lindquist, 1986; Huang et al., 
1999 ) , thus protecting cells against harmful insults. 
The genome sequencing of DaHSP23 in the present 
study showed that it is intronless, suggesting that it 
is inducible and qualified to be involved in diapause 
development and stress response. 

In D. antiqua, the expression of DaHSP23 was 
gradually up-regulated with SD and WD 
development, and declined to pre-diapause levels 
when SD and WD diapauses were terminated. The 
WD result was comparable to that found for HSP23 
in WD pupae of S. crassipalpis ( Yocum et al., 
1998). This is the first study on HSP23 expression 
in response to SD. Either cold- or heat-shock can 
upregulate the expression of inducible HSPs ( Yocum 
1998 ) ， 


evident for diapausing pupae that are not subject to 


et al., and up-regulated expression is 
temperature stress. This implies that expression of 
the sHSP protein is a component of the diapause 
syndrome in this species. Earlier studies suggest that 
small heat shock proteins are involved in the 
regulation of cell cycle arrest in diapause. Mounting 
evidence indicates that small heat shock proteins 
serve a pivotal role in cell growth and differentiation. 
If the small heat shock proteins are involved in cell 
cycle arrest associated with diapause, they are 
functioning in concert with other cell cycle regulator 
(Yocum et al., 1998). 

In contrast, the expression of HSP23 is not up- 
regulated as a function of diapause in D. triauraria 
and Lucilia sericata (Goto and Kimura, 2004; 
Tachibana et al., 2005). Thus, the up-regulation of 
the expression of HSP23 might not be a common 
feature of diapause insects. Variations among species 
are also found in the expression of other HSP genes 
such as HSP70 and HSP90. The expression of these 
genes was not regulated in D. triauraria and L. 
sericata during WD (Goto et al., 1998; Goto and 
Kimura, 2004; Tachibana et al., 2005), whereas 
that of HSP70 was highly up-regulated in D. antiqua 
during WD (Chen et al., 2006) and that of HSP90 
was down-regulated in S. crassipalpis (Yocum et al., 
1998; Rinehart and Denlinger, 2000). These 
results suggest that gene expression in response to 
diapauses might be species-specific. 
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A he FA a Be A DaHSP23 基因 的 克隆 及 在 
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摘要 :【 目的 】 低 分 子 量 (12 ~43 kDa) 热 激 蛋白 (sHSPs) 具有 抗 逆 应答 的 功能 , 滞 育 是 昆虫 抵抗 不 
良 环境 的 特殊 发 育 形式 ,但 sHSPs 在 昆虫 滞 育 发 育 过 程 中 的 作用 仍 不 清楚 。 本 研究 克隆 和 特征 化 
#88, Delia antiqua sHSP 基因 ,并 研究 它 在 夏 滞 育 和 人 冬 灌 育 发 育 过 程 中 的 表达 模式 ,为 阅 明 sHSPs 
在 滞 育 发 育 上 的 功能 黄 定 基础 。【 方 法 ] 通 过 RACE-PCR FARE T BH HSP23 基因 ,通过 相似 性 
比较 分 析 了 其 特征 、 结 构 域 及 与 双 翅 目 代 表 性 同 源 基因 的 系统 发 育 关 系 ;采用 实时 荧光 定量 PCR 
研究 了 该 基因 在 蔚 蝇 冬 滞 育 肾 和 夏 滞 育 肾 发 育 过 程 中 的 表达 情况 ,通过 表达 的 差异 比较 揭示 了 该 
基因 与 滞 育 发 育 的 关系 。 (AR) TERT BH HSP23 基因 ,命名 为 DaHSP23 (GenBank 登录 号 : 
HQ392521.1) ,其 cDNA 全 长 序列 为 904 bp ,编码 186 HAAR ,推测 蛋白 分 子 量 为 20.9 kDa, 等 电 
点 为 6.42。 该 基因 的 编码 蛋白 与 其 他 双 翅 目 昆 求 的 sHSPs 有 超过 66% 的 氨基酸 序列 一 致 性 ,与 已 
报道 的 其 他 双 翅 目 昆虫 的 滞 育 相关 HSP23 基因 同 源 。 基 因 组 测序 显示 该 基因 无 内 含 子 。 
DaHSP23 基因 在 黄 蝎 非 沿 育 蚂 的 发 育 过 程 中 一 直 保 桂 在 较 低 的 水 平 , 各 发 育 阶 段 间 的 表达 量 不 存 
JEL. LEA A Fe Lit ASAP KR Bai FACE AP ik BAS HARA, Bla Ra 
期 的 中 后 期 达到 峰值 ,在 滞 育 终止 期 逐渐 降 到 较 低 的 水 平 。【 结论 ] DaHSP23 AA EB MA iH A Fe 
夏 灌 育 发 育 过 程 中 明显 上 调 表 达 , 但 存在 差异 , 它 在 滞 育 期 的 调控 可 能 是 种 专 化 的 。DaHSP23 可 
能 在 昔 蝎 两 种 类 型 的 滞 育 上 起 重要 作用 。 

关键 词 : BH; HA; 热 激 蛋白 ; 基因 克隆 ; 表达 模式 ; 实时 荧光 定量 PCR 
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